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ABSTRACT 

An analysis of the steady and unsteady aeTodynamics of the space shuttle orbiter 
has been performed., It is shown that slender wing theory can be modified to account 
for the effect of Mach number and leading edge roundness on both attached and sepa- 
rated flow loads. The orbiter unsteady aerodynamics can be computed by defining 
two ecpjivalent slender wings , one for attached flow loads and another for th? vortex- 
induced loads. It is fotind that the orbiter is in the transonic speed region, subject to 
vortex-shock-boundary layer interactions that cause highly nonlinear or discontinuous 
load changes which can endanger the structural integrity of the orbiter wing and possi- 
bly cause snap roll problems . It is presently impossible to simulate these interac- 
tions in a wind tunnel test even in the static case. Thus, a well planned combined 
analytic and experimen’'al approach is needed to solve the problem. 
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Section 1 
INTRODUCTION 

The aim of the present analysis is to provide the building blocks for a computer 
program that can predict the aeroelastle characteristics of the space shuttle launch 
configuration through analytic extrapolation from es^erimental data, similarly to what 
was done for the ApoUo-Satum launch vehicle (Refs. 1-3). The most significant dif- 
ference between the present launch vehicle and the Saturn booster is the presence of a 
large lifting surfe.ce, the orbiter wing. As it is basically a delta wing, analytic means 
are needed for the prediction of the unsteady delta wing aerodynamics, including the 
effects of the leading edge vortices. Because of the asymmetric mating of the orbiter 
to the booster roll-pitch-yaw coupling effects become of concern and it is essential 
that the yaw-induced rolling moment on the orbiter delta wing can be predicted. 

Ihe first step towards the development of this needed analytic capability was 
taken in References 4, 5 and 6. In the present report this analytic development is 
taken one step further by including the vortex entrainment effect suggested in Refer- 
ence 5 and biy extending the analysis to account for Mach number effecte along the 
lines discussed in Reference 4. Finally, a critical look is taken at the orbiter wing 
aeroelastic stability charaeteristies without the complications introduced by the inter- 
ference flow field from HO Tanlc and SRM rockets . 
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Section 2 
DISCUSSION 

The anaiyees in 3Eleferences 4 and 5 will be repeated here to the extent necessary 
for good readability. The intent is to present the equations in a compressed, unified 
form such that the design engineer can get all heeded information frorn one source: 
this report. 

2 . 1 Delta Wing Aerodynamics 

hi incompressible flow, the nor jatal force on a slender, sharp edged delta wing 
can be determined as follows (Kef . 7): 

= Kp sin ca cos Gi + ICy. sin CH (1) 

Kp and Ky are parameters determining the magnitudes of attached flow and vor 

tox force components respectively, Ky is for all practical purposes a true constant, 

Ky « 7T, (Ref. 7); whereas Kp is more or less linearly dependent upon aspect ratio 

(Fig. 1), Jones slender wing theory (Ref . 8) applies only for = 1. The deviaticn 

between it and potential flow results (Ref , 7) at = 0 are represented as follows. 

It is assumed that the area denoted App in the inset in Figure 1 carries no load . 

This accounts in a crude manner for the trailing edge condition at M = 0 (as com- 

2 

pared to =1). It reduces the delta wing normal force by the factor cos 
.gbdng.. 

= !r(A/2)/ [l+ (A/#] (2) 
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This is the "First Approximation" shown in Fig. 1, The strip loading normal to 
the leading edge is 


(1/2) (a c„ /d^)_ = 




2 

7T sin 2 a sin 0- 


j^^/(b/2 c^) I 


(3) 


It was shown in Reference 5 that the potential flow loading on a delta wing never 
exceeded 75% of the maximum given by slender wing theory. Applying this "ceiling" 
changes Eq. (3) to 


(l/2Kd = 

Integrating Eq. (4) gives 


Tr ain 2 a sin 0^^^, 

(b/2o^) ~ * 


^ : I ^ 0.7 
0.7 > 0.7 


(4) 


= 0.91 sin a cos a. 
'ma " ^Na ^^a *' 

^ 0.04(1 -A|^TE> 


^^aTE 


^a ®LE 


(5) 


(77 = 4/3 TT for elliptie loading) 

Kp^ is given by Eq. (2). The value 0.91 is the "Second Approximation" in 
Fig. 1 which showa good agreement with potential theory for aspect ratios up to A = 3. 


It was also shown in Reference 5 that the vortex induced load distribution on a 
sharp-edged delta wing does not have the triangular shape prescribed by the conic flow 
assumption used in most theories, but has a "ceiling" similar to that for attached flow. 
The resulting load distribution for the 70% of the vortex induced loading that is located 
at the leading edge was shown to be: 


1 

1 . 72V d| A 


2 

( 1,72 7T ^ sin a: 

' 2 ' 

-<0*685 ;rsm a: 

...„2 


I '^0^4 
0.4 < ^ <0.9 
6.85 TTsin^ a (1-^) : 0.9 < ^ <1.0 


( 6 ) 
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and the rejmining 30% vortex induced loading produced by the vortex-entrainment 
effect over the center wing was shown to have the attached flow type distoibution given 
below 

1 = m 

0.3 2 \ / |o,735 ?r sin OJ; 0.7 < 1.0 

Integration of E^b. (6) and (7) gives: 

2 

Cj^ ^ 7T Sin £H 

^mV " ■'^^o^V^NV 

^^a ■' ^CG^] 

— 0.56 (1 - 7?t^ sin 0 j^-pi ) 

- 0.64 (l--^^sin^ 

= 4/3 7 T and rj^ is given by experimental data 
(see Fig. 2 and Ref. 9) 



It was shown in Reference 5 how the unsteady aerodynamics for the attached flow 
loads could be determined by using slender wing fiieory (Ref, 8) on an equivalent delta 
wing which has been shortened to provide the correct normal force defined by Fq. (5). 
The effective center chord for this equivalent wing is defined as follows 



0.955 cos 0 


LE 


—2 

2 - cos a 


11/2 



( 9 ) 
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Thus, the attached flow unsteady aerodynamles for oscillations in pitch (0) around 
^CG are given by slender wing theory in the following form 



From Fqs. (9) and (10) the center of pressure is obtained as 

I fsIL = 0.64(1-1 sin^ +—) (11) 

o 

whereas Eq, (S) gives 




= 0.64 (1 - 0.425 sin^* 


( 12 ) 


That is, the center of pressure is the same for slender wings, sin 0^^^ « 1, 
as it should be . 

Using Lamborae*s results (Ref. 10) it was shown in Reference 5 that the vortex 
induced unsteady aerodynamics could be determined as follows 


= -r I“-3 «Cg) + (iv- W] 

“v 


m 




■ -W 




0.3 sec Q!, 


V •- 




U 


~ ^CG xj ■ ^CG^ 


( 13 ) 
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wher& 1^/0 - 0*75 for oseiilailQns in pitch, and is given by Eq, (8); 

Fig. 3 shows the predictions by Eqs. (8) - (13) to agree well with experimental 
data (Refs, 11 and 12). 


2,2 Effect of Leading Edge Roundness 

Gersten has shown (Ref, 13) that leading edge roundness has a large effect on 
delta wing aerodynamics (Mg, 4), A large part of the measured decrease in lift and 
pitching moment is probably due to the delay of leading edge separation caused by 
leading edge roundness. Ihe 12.5% truncation of the wing tip (see shetch in Mg. 4) 
would not have any effect on the delta wing lift and moment according to the present 
analytic model (see inset in Fig, 1). Recent results by LaMar (Ref. 14) also indicate 
that the effect would be negligibly small. 


According to Ville (Ref. 15), leading edge separation should, in two-dimensional 




It is clear that the three- 


flow, occur at an angle of attack of o! 2 J }3 ~ \ o 

dimensional flow will delay the separation in the cross flow plane to an angle of attack 


“Ns > *^203 
ratio A - 

that C, 

Lmax 

obtains 


An estimate of this delay can be made by defining an effective aspect 
4 c^/b for half* the delta wing in the cross flow plane. It can be assumed 
is relatively independent of aspect ratio, (Refe, 16 and 17), Thus, one 


“Ns ^2Ds ^ ^LE^ 

For the small angles of attack of interest the; angle of attack o!_ is simply 
“s = “2Ds ^ ®LE - : 


*The leading edge flow is insensitive to what happens on the other wing half . 
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Experimental results (Eef^ 18) indicate that leading edge stall on the NACA-0012 
airfoil will occur in the angle of attack range 12.5“ < *^2Ds ~ depending upon 
the Reynolds number. For the orbiter main wing A = 45°, which in Eq. (15) gives 
ojg = o'gDg Thus, for =0 one obtains 17,7° < < 25,5“. 


The Mach number normal to the leading edge is simply 

ll/2 


,._2 


2 

cos Cl tan a + sin" 0^^ 


( 16 ) 


ITsing the two-dimensional data for the separated flow boundaries of the NACA 
64A 012 airfoil (See Ref. 19 and Fig. 5a) as a guide line, the prediction for = 0 
can be extended to transonic speeds as is shown in Figure 5b, The bars in Fi-^ure 5b 
represent oil flow data for the orbiter (Ref. 20). As the pictures were only taken at 
every 5 degrees of angle of attack, the bars span over Aa = 5“ , with the bottom in- 
dicating attached leading edge flow and the top of the bar established leading edge 
separation. Considering all the complications due to Qiree-dimensional flow effects, 
which will be discussed later in this report and in more detail in Part n (Ref. 21), the 
agreement is rather good between predicted and measured separated flow boundaries. 
It should be noted that without the aspect ratio correction, Eq. (14), the predicted 

Q! -values would be exactly half of what is shown in Figure 5b. For this reason the 
s 

predicted boundary in Figure 18 of Ref. 22 is much too low. 


For the test data in Figure 4 a value of suggested giving = 

3.95“. Substituting a with (o! - a ) in Eqs. (5) and (8) gives the rounded leading 
edge effect shown by the difference between the dashed and solid line curves in Fig- 
ure 4, Although the predicted effect of leading edge roundness is in the right direc- 
tion, it is less than what was observed experimentally. One possible reason for this 
is the following: It was discussed in Reference 5 how a slackening of the vortex-induced 
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load buildup was observed to occur aft of 40% center cbord (| > 0,4). This was 
attributed to a "loosening” of the vortex shedding; i.e. , the vortex becomes less con- 
centrated. Such an increase of the vortex core could occur for two possible reasons: 

1) the center core axial velocity is decreasing or 2) the vortex shedding mechanism 
from the leading edge has changed. If the second. reason is the significant one» the 
leading edge roundness could be expected to contribute to further "loosening” of the 
vortex shedding with associated loss in vortex-induced lift. However, if the ftrst 
reason dominates, which would be in accordance with the causative mechanism for 
the more severe ’'loosening" phenomenon called vortex burst (Refs. 23 - 25), the 
leading edge roundness effect should be accounted for by the (oi - o! ) correefion. 

This appears to be the case judging by the good agreement between predictions and 
experimental dynamic results (Ref , 11) for a delta wing with rounded leading edge 
(Fig. 6). 

That this purely static eftect of leading edge roundness suftices to explain the ■ 

dynamic effects is somewhat surprising in view of the large overshoot of static stall 
observed in dynamic testing of airfoils (Ref. 26). However, it is in agreement with | 

the dominance of three-dimensional flow observed by l^mboume (Ref, 10). He | 

showed that it is the flow conditions at the apex at an earlier time instant that deter- ; 

mine the instantaneous downstream vortex strength. 

2,3 Lateral Characteristics 

The dominant lateral characteristic of a Sideslipping delta wing is the roHing ^ 

moment derivative (Fig. 7). At an an^e of sideslip, the effective apex angle 
of the windward side is increased by an amount A 0^^ (see Fig. 8). i] 

-1 / i 

A^lj; “ ^ (tan IB/Gos a) (17) ^ 

or for small sideslip angles, ^ 

i 

. li 

A^Lg = j3 sec a (18) ij 

JL j 

For the 6% thick wing =12® gives =5.7° . ' 

. 2-7 , 'I 
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In attached flow, the windward side normal force is increasing because of the 
increased aspect ratio, Eqs. (1) and (2) give 


°N - ^ 

Kp = ZTTtaH 


( 19 ) 


Thus , the j3 -derivative for the normal force of the windward half of the delta 
wing can be determined as follows; 


1 1 . 


11 


2 C 




1%. 


de. 


LE 


d^ 


1 




= cot 9^^ - tan 


( 20 ) 


d0 


LE 


d|3 


sec a 


The corresponding rolling moment derivative for the wing half is 


1 ^ 

2 d|3 


a 

b 


2 dS b 




( 21 ) 


and the derivative for the full delta wing becomes 


•^^ta ^a’^a °Ha 

For the attached flow loads , the load distribution remains the same for small 
sideslip angles. This is not true for the vortex-induced loads, which axe determined 
only by the conditions at the leading edge. The loads given by Eq. (6) are the same 
regardless of whether it is a swept wing or a delta wing. That is, the presence of 
the center wing area of the delia wing has no efCeet on the vortex strength and the 
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vortex-iuduced loads near the leading edge^ given hy Eq.. (6), It is only needed to 
realize the vortex entrainment effect, Eq, (7). Thus, the windward half of the delta 
wing has the following increased load due to the sideslip angle (see Eig, 8). 


i = I (“■ “V® <“■ PVs C“. 0) 

s (0!, j3)/S (0;, 0) ~ 1 -i- /3 see a (cot 


(23) 


That is 


1 

2 V I 


2 cos a: 


(cot 0L^ -ian 


LE' 


(24) 


Eor the 70% of the vortex induced load located near the leading edge, Eq. (6), the 
effective lever arm jy for the rolling moment Is (see Fig, 8) . 


yy - ^y j^dy (tan ^ sec o;) 

- /3 sec O'! cos (jS sec o') 


(25) 


For small sideslip angles ^ the effective dimensionless lever arm for the wind- 
ward half of the delta wing is 


V 

1 - ??y ^COt 

A? 2 

^ Wy cos O! 


(26) 


For the 30% of the vortex load caused hy entrainment effects the dimensionless 
lever arm is ‘ 


111] = kJh 

\ b I.S 2 


(27) 
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The j3-derivative of the roUmg moment for the -windwarcJ half of the delta wing is 

1 _ 0.7 , 0.7 ^ a 

“2 d/3 " 2 d/3 \ b /.7 3 d§ \h /.7 

2 d/3 \ b /.3 


(28) 


Eqs, (24) - (28) define the following rolling moment derivative for the vortex in- 
duced loads on a delta wing 






d/3 7 /3^ “ 2 cos o! ^LE 


0.7 I 


V L 


(2 -tan 0^^)Tj^-l 


■t 0.3 (1 - tan 


(29) 


a 


The roll stabiliiy derivative given by Eqs , (22) and (29) in combination with 
Eqs. (5) and (8) is compared with experimental data (Eef. 27) in Mg. 9. Ihe agree- 
ment is excellent, substantially better than for Polhamus* theory (Refs. 6 and 27) 


JL 

•TT 


2.4 Mach Number Effects 


At sonic speed Jones’ slender wing theory (Ref. 8) applies. Ihus, the whole wing 
is effective, i.e,, = 0 and the ’’plateauihg" does not occur. The attached load 

distribution along the center chord is as follows for 


1 Na 


TT sln 2 a tan 6 


LE 


(b/2 c^) 


4 : 0 < ^ < 1.0 


2 d^ 

The corresponding vortex induced load distribution is 


(30) 


dC 

1 1 ^'"^NV 

0.7 2 df 


. 2 


1.72 TT I sin Cl ; 0 < f < 0.4 
0.685 7T sin^ Of : 0.4< ^ 1,0 


( 31 ) 


0^3 t “ 


0 < ^ < 1.0 




It should be noted that the present theory is of semi-empirical nature, whereas 
Polhamus’ theory is of die pure variety. 
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]iitegra.tion gives the following aerodynamic characteristics of a delta wing at 


= 1.0 

OQ 


= 7T (A/2) sin Q! cos a 


(«a - «Og) 


- - 2 2 „ 
^NV “ D! sec Q 


mv 


’a 


LE 

0.667 : 1^ = 0.587 


(32) 


An obvious way to make a smooth transition from = 0 to =1,0 in 

. . CO oo 

regard to is the following 


■^TE 


•TE‘ 


M^.=0 


M 


00 


(33) 


Applying the same smoothening also to the load distribution gives the following 
unified representation of the subsonic longitudinal aerodynamic characteristics of a 
delta'wing.'': 


% = 


= w (A/2) sin or cos a; ^1-0,09 V 1 - m/ ) 1 



/I * 2 

(1 - sm 

^LE 

4^ 


c = 

ma 

- (c./c) ( 



^CG> 


TTSitt^ Q! 

sec^ 

^LE 



/- ; 2 
11 - sm 

■^LE 

4i 


^mY "" 

- (P^/c) < 


‘^CG^ 


(34) 


(o. 


667 - 0.027 


•v/I-mJ) 


’’aS^^ ®IiE '^1"’'*^'/) 


“V 


(o. 

(l-)7 


Vl 


587 - 0.027 'VI -^ M 


■ 2 - 11 / 
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The unsteady aerodynamics are given by Eqs. (10) and (13) also at 5^ 0 if 
Eq, (9) is substituted by 

= 

Gq 

-sin" “V 1 - M^“) (2 - cos”" o!^^) J 
^N<rV’ ^a are obtained from Eq. (34) 

The expressions for the roiling moment derivative remain the same as before, 

Eqs (22) and (26). 

Figure 10 shows that the longitudinal aerodynamic characteristiGs of a delta wing 
at subsonic speeds predicted by Eq. (34) are in excellent agreement with expeinmental 
data (Ref, 27). Combining Eq. (34) with Eqs. (21) and (28) gives the predicted rolling 
moment derivatives at |3 ^ 0. Figure 11 shows that these are also in good agreement 
with Ihe experimental results (Ref. 27). 

It is of considerable interest to know how far into the supersonic region this 
approach based on Jones* slender wing theory (Ref. 8) can be carried. Brown (Ref. 28) 
Ms shown that the slender wing theory gives increasing over-prediction of the lift for 
increasing apex angle and increasing Mach number (Fig, 12). The restdts in Fig" 
ure 12 are well approximated by the following expression 




LOCKHEED MISSILES 8; SPACE COMPANY. INC. 



IiMSC-D0S7194 


K^ (m" ^ 1) ° ^ 


_ (^_^y 

\ tan |t / 


(37) 


When the leading edge goes supersonic, tan — tan |i, the vortex lift dis- 

appears, as it should. 

With these results the longitudinal aerodynamic characteristics of a delta wing can 
be deterndned as follows for supersonic speeds, > 1. 


Gw = 
Na 


ma 


'm 


tt Kj^ (A/2) sin a cos <r 

-(%/^H0-667- IcC3)Cna 


sin^ Q! sec^ ^ ^ 


ftan 9. 


tan 


(* / 


(38) 


•^mT = - <“o^) 


The unsteady aerodynamics are given by Eqs. (10) and (13) provided that Eq. (9) 
is substituted by 


o 


-2 , 
cos 


1 1/2 


(39) 


Combining Eqs. (37) and (39) with Eqs . (22) and (29) gives the rolling moment 
derivative for supersonic speeds. 

At high supersonic speeds , where the leading edge flow is supersonic, another 
analytic approach is needed. The unsteady embedded Newtonian flow theory (Eefs. 29 
and 30) can provide the basis for prediction of the attached flow loads, on top of which 
separated flow effects can be superimposed (Ref. 24) . 
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2.5 EJEfect of Trailing Edge Sweep 

The delta wing analysis can be extended to apply to slender wings with swept 
leading and trailing edges. Two equivalent delta wings are defined, one for the 
attached flow loads and another for the vortex indueed loading (see Fig. 13) . The 
equivalent delta wing characteristics are referenced to the slender wing geometry 
using the following ratios^ 

SVS = (1 “ tan e^^tan -7? tan 0^^ tan 

= 1/(1 - 77 tan tan 

= (1 - tan tan - rj tan tan 0^^ 

Figures 14 and 15 show that the low speed longitudinal and lateral stability char- _ 

acteristies predicted by use of Eq. (40) in combination with the delta wing expressions, i 

Eqs. (4), (7)i (22), and (29), agree well with experimental results (Hef. 27). It is 
apparent that the present method offers a definite improvement over the Polhamns ' [ 

theory, especially in regard to the arrow wing characteristics . This is true also at 

high speed, =0.8 (Figs. 16 and 17). r7j 

■ ■ ' i 


■ V '1 
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The star indicates parameter values for the equivalent delta wing, 
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Section 3 

OHBITER UNSTEADY AERODYNAMICS 

space shuttle orbiter has a wing of double-delta pianform (Eig, 18) , An 80“ 
swept inner delta wing, a strake, is preceding the 45“ swept main delta wing. It is 
well known that the strake vortex induces a load on the main wing (See Fig. 19 and 
Refs. 31-33) . At high angles of attack the strake and main wing vortices may combine 
to form one (large-core) vortex (Ref. 32). In absence of a fuselage the downstream 
sirength of the strake vortex, determining the ma^itude of tlie vortex induced loads, 
would be the integrated result of the angle of attack distribution along the strake lead- 
ing edge. However, in the presence of the orbiter fuselage (Fig. IS) the flow situa- 
tion changes dramaticaEy. At low angles of attack a comer separation occurs at the 
wing-fiiselage juncture which reaches all the way to the strake apex already at moder- 
ate angles of attack. The comer separation is vented via a vortex which affects the 
lift over the aft wing. When the separation oceucs at the strake apex, the situation is 
similar to that for ficee-body vortices on slender bodies of revolution (Refs, 34, 35, 

38 and 39). For the shuttle this means that crossflow at the strake apex determines 
the proximity of the vortex core to the aft wing, thereby determining the aft wing lift. 
Thus the situation is very much different from the one dealt with earlier in the pure 
delta wing analysis. 

3,1 Rigid Body Dynamics 

The rigid body dynamics of the orbiter (Fig, 18) are computed in the following 
manner. An equivalent slender wing is defined for computation of the attached flow 
unsteady aerodynamics, similarly to what was done in the delta wing analysis earlier. 
The trailing edge of the equivalent Wing is located such that the computed slender 
wing force derivative at cj = 0 agrees With the measured 
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induced loads are defined as the difference at angle of attack between actually measured 
static characteristics (iiefe, 36 and 37) and the computed attached flow characteristics. 
In this manner the flow complications caused by the fuselage are accounted for in re- 
gard to thp- magnitude of the vortex induced loads. 3h order to obtain the unsteady 
aerodynamics one has to determine the phasing of these loads. 

*From the earlier discussion it appears reasonable to lump the crossflow effects 
on the stralce-fuselage vortex to the strake apex. That is , the strake-fuselage vortex 
is lil?e a free body vortex (Refs. 34, 35, 38 and 39), or a partial span vortex (Refe, 4 
and 40), and one can assume that R = XJ^ . Thus, the lift generated by the straJce- 
fuselage vortex at a station x-x^ downstream of strake apex is determined by the 
cross flow at the strake apex a time instant At earlier, where At = (x-x.)/XJ . 

xi. CO 

This analysis is carried out in detail in Part It of the present report (Ref, 21). The 
dynarndLc characteristics computed in this manner are shown in Figures 20 a: id 21. 

The agreement with experimental results (Ref. 41) is good, not only at subsonic speeds 
(Fig, 20) but also in the transonic speed region (Hg. 21). Note the opposite effects on 
dynamic and static stability of the separation induced loads. In the transonic speed 
region (Fig, 21) shock-boundary layer interactions complicate the picture. However, 
also the shock-induced flow separation is to a large extent controlled by the strake- 
fuselage vortex (Ref, 21), Thus, the cross flow at the strake apex determines also 
the shock-induced separated flow effects on the vehicle dynamics . 

3.2 Aeroelastio CharacteristLes 

The rigid body dynamic data for the orbiter shown in Figure 22 for - 1.2 
reveal that a dramatic change of flow pattern occurs at cu « 8“, The measured large 
increase of dynamic stebnity and corresponding moderate decrease of static stability 
are typical for nonlinear, possibly discontinuous, aerodynamic characteristics which 
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are associated vidtli suddenly increased flotv* separation, such as has been observed on 
slender cone** cylinder bodies (Hef. 42). This acconated for the 1% (of critical) loss of 
damping measured on the Saturn- Jupiter nose shroud (See l?ig, 23 and Refs. 42 and 
43). Whether this sudden increase of the flow separation will increase or decrease 
elastic vehicle damping depends entirely on the mode shape. Figure 24a shows how one 
of the candidate straight wings for the early space shuttle configuration experienced 
xindamping (negative damping) at - 0.85 (Ref, 44). Figure 24b shows the meas- 
ured vibration response for a slightly different wing. However, on the swept wing 
of a high performance aircraft a sudden change of the shock-induced separation pat- 
tern causes increased damping (Fig. 25 and Ref. 46). Ihat the wing bending response 
still increases in this case is the result of the increased forcing function, the buffet 
input. It is the combined effect of the forcing function and the change of aerodynamic 
damping that determines the buffet x’esponse (Ref. 47) 


Figure 26 shows the similar effects of sudden separation on cone-cylinder bodies 
(Ref. 48) and two-dimensional airfoils (Ref. 19). A force couple is generated together 
with a net negative normal force, which may be more pronounced in the two-dimensional 
flow ease. Figure 27 illustrates how this net negative force would produce a statically 
destabilizing and hence damping effect on the rigid body orbiter, oscillating in pitch, 
all in agreement with the results shown in Figure 23. The effect for the elastic mode 
sketched in Figure 27 wiH be the opposite with the force couple providing a statically 
restoring and hence undamping moment. It should be noted that it is the force couple 
that does the damage to the aero elastic damping. Thus, if the resultant net force had 
been of negligible magnitude the rigid body dynamic data would not have ^ven any 
warning about this possibiliiy of aeroelastic instability. The question is then, whether 
or not the actual orbiter can experience such large adverse aeroelastic effects. Flow 


visualization photographs (Ref. 49) indicate that this can indeed be the case. The 
nodal line for the first torsional mode, as obtnined from Reference 50 , is delineated 


"^Actually the fin for the same straight wing space shuttle configuration, (Ref, 45). 
The pealc response is, therefore, at = 0.90 and not at = 0.85. 

^Note that the flutter investigation in Reference 50 was perfornied at ch = 0 and 
(X = 3“, i.e. , far from the critical a-range for sudden separation changes. 
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ia figure 28. One notices that the separation line moves from a position aft of the 
nodal line to a position &r forward of the nodal line when the angle of attack is in- 
creased from O' — 5° to o; = 10“. 


When combining the information presented in S'igures 22 throtigh-28 one becomes 
CGHvineed that large adverse aeroelastie effects will be experienced by the or biter 
wing at certain critical combmations of angle of attack and ,Mach num^ 
shows how this critical angle of attack varies with Mach number and cone angle for 
cone-cylinder bodies (Refs. 47 and 48). Within the oi^^^^-range shown the flow alter- 
nated between the two separated flow patterns . TEiis would indicate fliat the adverse 
aeroelasiic effect can be realized for a range of angles of attack a - * "^“orit * 

as is also verified by the results obtained with the Jupiter nose shroud (Mg. 23), 

^crxt ^ particular, are very sensitive to Reynolds number. One can, 

of course, expect the critical angle of attack region for the orbiter wing to be equally 
sensitive to Reynolds number * ha addition, eleven deflection will strongly influence 
oicj,it for the orbiter . 


The orbiter wing with its highly swept leading edges is subject to strong three- 
dimensional flow effects. How the wing sweep affects tiie shock-boundary layer in- 
teraction has been studied thoroughly (See Ref. 51). The similarily with the flow 
pattern on the orbiter wing (Ref . 49) is apparent in Fi^re 30. hi addition to this 
triple-shock complication, the wing sweep also introduces a unique boundary layer 
transition behavior with associated decisive effects on the flow separation patterns. 

The spanwise flow, with its propensity for the development of an inflexion point in its 
boundary layer profile (Ref. 52), causes the flow to beconie turbulent on the outer wing 
panel while it is still laminar inboard. For some ranges of Mach number and angle 
of attack this causes the outboard wing to have a region of attached leading edge flow 
with retarded shock-induced separation (see Mg, 31 Ref , 53). ^ At a certain angle 

of attack the separation jumps to the leading edge causing a dramatic change of the 
flow Held over the complete wing. How this critical angle of attack varies with, lead- 
ing edge sweep for TACT~F111 (Refs, 54 and 55) is shown in Mgure 32. The Fill 
wing with its supercritical airfoil and built-in 4 degrees ''wash-out" at the wing tips 
does experience the flow pattern shown in Mgure 31, whereas the orbiter with its 
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more conventional airfoil sliape does not, (Fig. 28). TliatiB, on. tlie or biter the 
separation does not jump all the way to the leading edge wben o; exceeds 
Othervrise, orbiter fits into the general trend shown by the TACT- Fill 

(see Fig. 32). 


Unfortunately there are no flow pictures available for the orbiter without 03VIS- 
pods . Without the restraining action from the pod shock the separation probably 
moves farther aft before it jumps forward. Whether or not it jumps all the way to the 
leading edge is a moot question. 3h any ease, the discontinuous force change will he 
larger than in the presence of the OMS-pods, This explains the much larger damping 
measured on the "clean’’ orbiter (see Kgi 22), It is also in agreement with the ob- 
seiwed large effect of pod geometry on the complete wing flow field (see Fig. 33 and 
Ref. 49). It is obvious that at a sideslip angle will be different for the up- 

wind and down-wind wing halves . That is, the snap roll problem encountered on the 
straight-winged early orbiter configuration (Refs. 56-58) could be of concern also for 
the present orbiter . The wing tip loads dominate the rolling moment and have an 
even larger effect on the wing root bending moment (Ref, 59). hi addition to sideslip, 
control deflections can have very large effects on elastic and rigid body dynamics 
(Refs, 27, 35 and 56). During launch the interference flow effects from HO-tank and 
SRM-rockets are large even when not considering the effects of rocket plumes (see 
Ref. 60 and Vol, IE of the present report). Adding to all this complexily is the fact 
that it is very difficult to simulate shock boundary layer interaclioh, even in the simple 
case of a straight wing (see Fig, 34 and Refs. 61-63), When considering the three- 
dimensional aspects (Figs. 31-34) the problem becomes even more formidable. 
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Section 4 
CONCLUSIONS 


A study of the steady and unsteady aerodynamics of the space shuttle orbiter has 
given the following resulte . 

o The effect of Mach number for a subsonic leading edge is accounted for by a 
simple modification of dbnes* slender wing theory for the attached flow loads 
and by a reformulation of Polhamus' theory for the vortex induced loads. 

9 The effect of leading edge roundness is largely accounted for by considering 
the delay of leading edge sepa.i:ation in the cross flow plane due to leading 
edge roundness. 

0 The effect of trailing edge sweep (forward or back) is well predicted by use 
of two equivalent delta wings, one for the attached flow loads and another for 
the vortex induced loads. 

• Both pitch and roll stability derivatives, determined by the presented closed 
form solutions, are in excellent agreement with available experimental data. 

• The slender wing analysis can be extended to the orbiter wing with its double- 
delta plan form wing by defining two equivalent wings. 

The attached flow loads are given by an equivalent slender wing that 
gives the measured C^^ at O! = 0. 

The vortex induced loads are defined as the difference at a > 0 between 
measured total loads and computed attached flow loads. The cross flow 
at the strake apex is controlling the unsteady effects of the vortex in- 
duced loads. 

• The unsteady aerodynamics of the orbiter computed in this manner are in 

good agreement with measured dynamic characteristics in the whole Mach 
number range investigated, 0.3 < “ 1*2, 
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« ]h iJie transoEic speed region, 0 . 9 :< 1,2, the space sliuttle orbiter 

is subject to yorteic-shock-boimdary layer interactions which cause highly 
nonlinear or discontinuous changes of the aerodynamic loads when a critical 
angle of attack is exceeded. It is found that this could have a strongly ad- 
verse effect on the aero elastic stability of the wing in its lowest torsional 
mode. This flow phenomenon can also have a dramatic effect on the vehicle 
dynamics: with the possibility of snap roll being of some concern. 

« The vortex-shock-boundary layer interaction is extremely sensitive to model 
surface roughness and free stream Beynolds number, and it appears impossi- 
Ide to simulate it in a wind tunnel test unless full scale Reynolds number can 
he reached. This is not possible with present ground test facilities, at least 
not with the size model needed for the complicated space shuttle, Por these 
reasons the following approach is suggested. 

Obtain consistent static and dynamic data for a "nominal current" con- 
figuration including "current" OMS-pods and the operational range of 
control deflections. 

Conduct parametric tests to determine incremental effects of geometric 
changes, control deflection, sideslip, etc, as well as of surface rough- 
ness and Reynolds number. 

Develop analytic means by which the wind tunnel test results can he 
predicted thereby ensuring confident analytic extrapolation to full scale 
conditions. 

Obviously the investigation outlined above should be broad enough to enable 
the analyst to find "fixes" as needed to ensure structural integrity and 
acceptable vehicle dynamics. 
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Appendix A 
NOMENCLATURE 


A 

a 

“^TE 

b 

c 

% 

Kp, Ky 
Kpi 

L 

I 

M 

M 

P 

N 

P 

q 

^N 

Be 

S 


2 

aspect ratio, A = b /S 
speed, of sound 

inefficient wing area at M = 0 (Fig. 1) 

wing span 

reference length 

slender wing root chord 

Mach number parameter, Eq. (36) 

potential flow and vortex lift :fe.ctors, Eq. (1) 

first approximation of Kp 

lift: coefficient = L/(p^ U^/2)S 

rolling moment; coefficient ufj/2)Sb 

Mach number: M = U/a 

pitching moment: coefficient C„ = M /(p U /2)Sc 
^ m p <ao 

normal force: coefficient ^ 

2 

static pressure: coefficient = (p - P^)/(P^ 

pitch rate 

airfoil nose radius 

Reynolds number 

reference area (= projected wing area) 
local semi-span 
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t 

U 

u 


X 


y 

a 


a 


o 




5 (x, t) 


V 

e 


LE 


TE 


I 

P 

O' 

0(X) 

CO 

CO 


time 

horizontal velocity 
convectioii velocity 

axial body-fixed coordinate (see inset in Fig. 1) 
spanwise body-fixed coordinate (see inset in Fig. 1) 
angle of attack (Fig. 6 ) 
trim angle of attack 
sideslip angle (Fig. 6 ) 

elastic vehicle deflection, 5 (x, t) - <^(x) q(t) 
damping, fraction of critical damping 
dimensionless y-coordinate, 77 = y/s 
angular perturbation in pitch 
cone half angle 

apex half angle (see inset in Fig. 1) 

trailing edge sweep angle (Fig, 12) 

leading edge sweep angle,, A= tt/2 - 

-1 

Mach angle, y. = cosec (M^) 
dimensionless x-coordinate, ^ = (x^ - x)/c^ 
air density 
wing tip acceleration 

x-distribution of normalized bending deflection, 5 (x, t) = <^(x) q(t) 
free-free bending frequency and rigid body pitching frequency 
reduced frequency, co = co 
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Subscripts 
A. apex 

a attached flow 

CG center of gravity 

crit critical 

d downstream 

eff. effective 

LE leading edge 

max maximuni 

N or X normal to leading edge 
s separated flow 

TE trailing edge 

V vortex 

2D two-dimensional flow 

CO freestream conditions 

Superscripts 

(*) trailing edge coordinate, Eq. (39). 

(-) barred quantities denote integrated mean values, e.g. , centroid of 

aerodynamic loads 
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Figure 1 Attached Flow lift Factor of Delta Wings at = 0 
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I’igure 4 Comparison of Predicted and Measured Effects pf Leading 
Edge Roundness on Static Characteristics of a Slender 
Wing of Aspect Ratio A = 0.78 
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Figure 5 Flow 
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bo Flow Separation Boundary of the Space 
Shuttle Orbiter 
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Figure 9 Comparison Between Predicted and Measured 
Yaw Stability Derivative 
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Figure 10 Effect of Mach Number on Longitudinal Delta Wing Aerodynamics 
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Figure 13 Definition of Trailing Edge Parameters 
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Figure 20 Orbiter Dynamics at Subsonic Speeds 
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a. Measured Damping 

Figure 24 Bending Response and Damping Characteristics of 
Straight Wings at High Subsonic Mach Numbers 

F-24 


LOCKHEED MISSILES 8t SPACE COMPANY. INC. 


LMSC-D057194 



INTENSITY WITH 



b. APPARENT DAMPING OF 1ST 
mNG-BENDING MODE 


Figure 25 Bending Response and Apparent Aerodynamic Damping of the 
Swept Wing of a Hi^ Performance Aircraft at 0.7 
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Figure 27 Sketch of Similar Shock-Induced Force-Couples on 
Cone- Cylinders and Airfoils 
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Figure 29 Critical Angle-of-Attaclc Regions for Cone-Qylinder Bodies 
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Figure 33 Effect of OMS-Pod Geometry on Orbiter Wing Flow 
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Figure 34 Scale Effects on Shock- Boundary Layer Interaction 

at M =0.8b 

CO 
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